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RATE OF CONVERGENCE IN THE MULTIDIMENSIONAL 
CENTRAL LIMIT THEOREM FOR STATIONARY PROCESSES. 
APPLICATION TO THE KNUDSEN GAS AND 
TO THE SINAI BILLIARD 

By FRANgoiSE Pene 
Universite de Bretagne Occidentale 

We show how Rio's method [Probab. Theory Related Fields 104 
(1996) 255-282] can be adapted to estabhsh a rate of convergence in 

in the multidimensional central limit theorem for some stationary 
processes in the sense of the Kantorovich metric. We give two appli- 
cations of this general result: in the case of the Knudsen gas and in 
the case of the Sinai billiard. 

0. Introduction. 

0.1. Context. We denote probability dynamical system (fi, J^, T) where 
z^) is a probability space endowed with a ;/-preserving transformation 
T of $7. Let a probability dynamical system z^, T) and a measurable 

function f ■.Q.^W^ (with d > 1) be given; we consider the stationary process 

We say that a M'^-random variable N is Gaussian if, for any /3 G M*^, the 
distribution of the real-valued random variable {(5, N) is either a normal dis- 
tribution or a Dirac measure. With this definition, a Gaussian random vari- 
able has a general normal distribution (cf. [16], III-6 for more details). We say 
that we have a central limit theorem (CLT) for (Xfc)fc>i if (-^ Z]fc=i ^fe)n>i 
converges in distribution to a Gaussian random variable A^. 

CLTs in the context of dynamical systems have been established in many 
articles (cf. [13, 27, 36, 39]). In these works, multidimensional central limit 
theorem follows directly from one-dimensional central limit theorem. In- 
deed, let us recall that the fact that (-^ Z)fc=i ^fc)n>i converges in dis- 
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tribution to N means that, for all (3 £W^, the one-dimensional process 
(:^Efc=i(/?>^fe))n>i converges in distribution to {(3,N). 

In the present paper, we are interested in questions of speed of convergence 
in the CLT for multidimensional stationary processes. There are many ways 
of estimating the speed in the CLT. Let us endow M*^ with the supremum 
norm | • |oo and with the Borel cr-algebra B{M.'^). For any (p:W^ ^ M, we 
denote by its Lipschitz constant: 

:= sup — , _ , . 

We can estimate the rate of convergence in the CLT for {Xk)k by estimating 
the following quantities: 

(a) uniform norm of the difference between the distribution functions 
{DF metric): 



DFn ■■= sup 

{xi,...,xa)& 



-in n 



k=l 



k=l 



where xj^^ and N^^^ are the ith coordinates of X^ and of N, respectively; 
(b) in the sense of the Prokhorov metric (H metric): 



n„ :=inf<^ e>0:VSG 



— J2Xk€B] - F{N€B')< 



k=l 



where is the open e-neighborliood of B; 

(c) in the sense of the Lipschitz bounded metric {LB metric): 



LBn :=sup-^ 



1 



fc=i 



■E[cf>{N)] 



(d) in the sense of the Kantorovich metric (k metric): 



Kn :=SUp-^ 



k=l 



mm 



We will give additional details about the metrics corresponding to these 
quantities. 



0.2. Previous results. 
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0.2.1. One- dimensional processes (d = 1). When d = 1, it is classical to 
estimate the speed in the CLT in the sense of the uniform error between the 
distribution functions {DF metric). 

In [1, 3, 15] a rate of convergence in in the sense of the DF met- 
ric is established for sequences {Xk)k of independent identically distributed 
random variables such that E[|Xip] < +oo. Moreover, this result is optimal. 

This result has been extended to some martingale processes (cf. [7]) and 
to some stationary processes ([37] extended in [22, 23] and in [25, 26]). 

When d = 1 and when (Xfc)fc is a sequence of independent identically 
distributed random variables, Nagaev [31] establishes a nonuniform estimate 
for the difference between the distribution functions: 



A direct consequence of this is a speed of convergence in in the sense of 
the DF metric but also in the sense of the Kantorovich metric (cf. Proposi- 
tion 0.10). Results in the sense of the Kantorovich metric have been estab- 
lished by many authors. Let us mention the article [41] of Sunklodas. 

0.2.2. Multidimensional processes. For sequences of independent identi- 
cally distributed random variables (^fc)fc>i with values in and admitting 
an invertible covariance matrix and admitting moments of the third order, 
a speed in is established by Bergstrom in the sense of the DF met- 
ric (cf. theorem of page 121 in [2]). This result gives an extension of the 
Berry-Esseen result to the d-dimensional case. 

In [23] Jan shows that Rio's result can be extended to the d-dimensional 
case (in the sense of the DF metric). 

In [44] Yurinskii establishes an inequality linking the Prokhorov metric 
with characteristic functions. This result allows to establish a rate of con- 
vergence in -A= for sequences of independent identically distributed random 
variables in the sense of the Prokhorov metric. We will recall and use this 
inequality in the case of the Knudsen gas. 

Let us also mention the works of [4, 35, 38, 42] in which the rate of 
convergence in the CLT is estimated in other ways for independent random 
variables sequences. 

0.3. Contents of the present paper. In Section 1 we give a result of speed 
of convergence in the CLT for some stationary processes {Xk)k in the sense 
of the Kantorovich metric (Theorem 1.1). The proof of this result, given 
in the Appendix, uses an adaptation of the method developed in [37] and 
extended in [25, 26]. In these papers, a rate of convergence in the sense of the 
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DF metric has been established in the one-dimensional case. Our hypothesis 
is analogous to the hypothesis of [25, 26] but weaker than it. 

In Sections 2 and 3 we will give applications of our result. 

In Section 2 we study the Knudsen gas model studied by Boatto and 
Golse in [6]. We use a Markov model of it. We show that the results of [6] 
are related to questions of rate of convergence in the CLT. For this model, 
we will, first, estimate the speed of convergence in the CLT in the sense of 
the Prokhorov metric (using Yurinskii's result of [44] and an extension of 
the method of perturbation of quasi-compact operators [20, 29, 30]). Second, 
we apply Theorem 1.1 of Section 1 and establish a rate in the sense of the 
Kantorovich metric. This result gives an extension of a result of [6]. 

In Section 3 we are interested in the question of the rate of convergence in 
the multidimensional CLT for {Xj. := f oT^)^ where T is the billiard trans- 
formation of the Sinai billiard [40] and / is a smooth (Holder continuous) 
function. Using Theorem 1.1 of Section 1, we establish a rate of convergence 
in A= in the sense of the Kantorovich metric. This is, to our knowledge, the 

first time that a speed of convergence in is established in this context. In 

[33, 34] a speed of convergence in ^^/^-q for any a > has been established 
in the sense of the DF metric and in the sense of the Prokhorov metric 
(with an adaptation of a method developed by Jan in [23] using characteris- 
tic functions). The result we present here does not exactly improve [33, 34]. 
It gives a better rate with another metric. 

0.4. Some metrics for probability measures on . Let us denote by 
A^i(R'^) the set of probability measures on (R'^', ^(M'^)), where B{W^) is 
the Borel fi-algebra of W^. 

Definition 0.1 (The DF metric). For all P,Q in MiiW^), we define 



Definition 0.2 (The Prokhorov metric; cf. [5, 14]). For all P,Q in 
Mi{W^), we define 

n(P, Q) := inf{e > O:\fBe B{R'^), (P(S) - Q(-B^)) < e}. 

Definition 0.3 (The Ky Fan metric for random variables). If X and Y 
are two M'^-valued random variables defined on the same probability space 
(£'o,'?OiIPo)i we define 




JC{X, Y) := inf{e > 0:¥o{\X -Y\^ > e) < e}. 
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Let us recall that linin^^oo K,{Xn,Y) = means that {Xn)n converges in 
probability to Y. 

Proposition 0.4 ([14], Corollary 11.6.4). For all P,Q in Mi{W^), the 
quantity Il{P, Q) is the infimum of}C{X, Y) where X and Y are two Mf^ -valued 
random variables defined on the same probability space and such that the dis- 
tribution of X is P and the distribution ofY is Q. 

Definition 0.5 (The bounded Lipschitz metric). For all P, Q in A^i(M'^), 
we define 

BL{P, Q) := sup{|Ep[0] - Eq[,/)]U : M"^ ^ M, +L^< 1}. 

In particular, for any bounded Lipschitz continuous function ^M, 
we have 

\Ep[4>]-Eq[4>]\<BL{P,Q) X (11011^ +L^). 



Proposition 0.6 ([14], Theorem 11.3.3). Let {Pn)n he a sequence of 
AiiiW^) and let P be in AiiiW^). The following properties are equivalent: 

(i) the sequence {Pn)n of probability measures converges weakly to P; 

(ii) lim„^+oon(P„,P) = 0; 

(iii) lim„.^_|_oo 

(iv) lim„,_>+oo ^^(-fn) -P) = 0; if P has a continuous distribution. 

More precisely, we have (cf. [28], Proposition 1.2 and [14], Problem 11.3.5): 

lBL{P, Q) < U{P, Q) < ClBLiP, Q)f\ 

Let us denote by Ali^int(K'^) the set of probability measures on (M*^, B{W^)) 
admitting moments of the first order. 

Definition 0.7 (The Kantorovich metric, cf. [14, 16]). For all P,Q 
m 7Wi,int(M'^), we define 

k{P,Q) :=sup{lEp[(/.] -Eq [</>]], (/.iM'^^M, < 1}. 

In particular, for any Lipschitz continuous function (^irM*^ — > M, we have 
|Ep[,^]-Eq[(/)]1 <k(P,Q) xL^. 
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Proposition 0.8 ([14], Theorem 11.8.2). For all P,Q in MiiR"^), the 
quantity k{P,Q) is the infimum o/E[|X — l^|oo]; where X and Y are two 

-valued random variables defined on the same probability space and such 
that the distribution of X is P and the distribution of Y is Q. 

Proposition 0.9. Let {Pn)n be a sequence of Mi^mti'^''') and P in 
-A4i,int(IK'^)- The following properties are equivalent: 

(i) the sequence {Pn)n of probability measures converges weakly to P 
and we have lim„^+oo /^d \x\oodPn{x) = /jjd \x\oodP{x); 

(ii) lim„^+ooK(^'n,^') = 0. 

Proof. According to Proposition 0.6 and to the fact that BL{Pn,P) < 
K{Pn,P), it is easy to see that (ii) impHes (i). 

Let us now suppose that (i) is true and prove that (ii) is then true. Let 
us write a„ := \ J^d \x\oodPn{x) — J^d \x\oodP{x)\. Let (/):]R'^ — > M be any 
Lipschitz continuous function with Lipschitz constant bounded by 1. For 
any nonnegative real number M, we define : I^'^ — ^ by 

(<p{x), if\(t>{x)-m\<M, 

^m{x) = { ^{<d) + M, if (l}{x) > (/.(O) + M, 
[ 4){<d) - M, if <j){x) < (/.(O) - M. 

For ah M > and ah integer n > 0, we have 

\EpJ^PM]-^p[^pM]\ < (M + 1)5L(P„,P) 

and 

VX G M"^ \i^M{x) - 0(x)| = I^m(x) - (A(x)|l{|^|^>M} 

< (|x|oo -M)1{|^|^>M} 

and therefore 

IEpJVm -^]\< EpJd • loo - M)1||.|^>M}] 
<EpJ|-|oo-min(|-|oo,M)] 

<Ep[\ ■ \^] + an - Ep[mm{\ ■ \^,M)] + {M + l)BLiPn,P) 

< Ep[(| • U - M)1||.|^>M}] + an + (M + l)BL{Pn,P). 

Hence, for any M > and any integer n > 0, we have 

|Ep„ [0] - Ep[0]| < 2Ep[(| • loo - M)1||.|^>M}] + a„ + 2(M + 1)5L(P„,P). 

Let a real number e > be given. Let us fix such that Ep[| • |ool{|-|oo>J\/e}] ^ 
|. Let be such that, for any integer n > N^, we have ^ f and (Af^ + 
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l)BL{Pm P) < I (according to Proposition 0.6, such a A^^ exists). There- 
fore, for any n> N;, and any Lipschitz continuous function (/> : M'^ ^ M with 
Lipschitz constant less than 1, we have 

|Epj0]-Ep[(/.]|<e. □ 

Proposition 0.10 (The Kantorovich metric, case d=l ([14], problem 
11.8.1)). For all P, Q in A^i,int(M), we have 

k(P,Q)= / \P{]-oo;x]) - QQ-oo;x])\dx. 

JR 

1. Abstract theorem. We will denote by Lip(M'^,M) the set of Lipschitz 
continuous functions from M"' into M. Let a probability space {Q,,J-,u) be 
given. For any i/-integrable function /lO— >M, we denote by ^u[f] the ex- 
pectation of / with respect to probability measure i^: 

E,[/] := / fdu. 
Jn 

For all real- valued functions f,g in L^(0,z/), we recall the definition of the 
covariance of / and g (with respect to measure v): 

CoyM, g) = ^u[fg] - ^uimAg]. 

For all A = {ai, . . . , ad) and B = {bi, . . . ,bd) in M*^, we denote hy A (S) B 
and j4®^ the square matrices given by 

Ai^ B := {aibj)ij=i^,„^d and A'^'^:=A^A. 

Let M = (-Mjj)jj=i^...^rf be a random variable on with values in the set 
of the square matrices such that, for any i,j = l,...,d, Mij is z/-integrable. 
Then, the expectation E,y[M] of M is the d-dimensional matrix given by 

E,[M]:=(E,[M„-])ij=i,„„rf. 

Let us consider a sequence of stationary M'^-valued random variables 
{Xk)k>o defined on (0,JP", i/). 

For any integer n > 1, we write 5„ := J2k=i -^k and Sq = 0. 

Using an adaptation of Rio's method developed in [37], we will establish 
the following result: 

Theorem 1.1. Let (Xfc)A:>o o, sequence of stationary W^-valued 
hounded random variables defined on {Q,J-',u) with expectation 0. Let us 
suppose that there exist two real numbers C >1, M > max(l, ||Xo||oo) o.nd 
an integer r > and a sequence of real numbers {<fp^i)p^i bounded by 1 with 
X]p>iPniax;=o,...,[p/(r+i)J V^p,/ < +c« such that for any integers a,b,c>0 
satisfying l<a + b + c<3, for any integers i,j, k,p, q, I with I < i < j < 
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k<k+p<k + p + q<k+p + l, for any ii, ^2, ^3 £ {1, • • • > d}, for any 
F -.W^ X ([— M;M]'^)^ bounded, differentiahle, with hounded differential, 
we have 

<C(||F||l^ + \\\DF\oo\\L^)iPp,u 

is) 

where DF is the Jacobian matrix of F and Xm is the sth coordinate of Xm- 
Then, the following limit exists: 

T?:= lim -nsT]). 

IfY? = d, then the sequence {Sn)n is bounded in . 

Otherwise the sequence of random variables i^)n>i converges in distri- 
bution to a Gaussian random variable N with expectation and with co- 
variance matrix and there exists a real number B > such that, for any 
integer n > 1 and any Lipschitz continuous function ^ M, we have 



E 



■E[0(iV)] 



n 



The proof of this theorem is given in the Appendix. 

2. Application to the Knudsen gas. Fohowing Boatto and Golse [6], we 
are interested in a generahzed model of the Knudsen gas with an isotropic 
component. In the present section we use a probabilistic approach. We show 
how this problem can be modeled by a Markov chain. Using the method 
of perturbation of operators due to Nagaev (see [19, 20, 29, 30]), we get 
a rate in in the multidimensional CLT in the sense of the Prokhorov 

V". 

metric. Moreover, we establish the same rate for the Kantorovich metric. 
This second result is an application of Theorem 1.1 and gives an extension 
of a theorem of [6]. 

2.1. The model. In this section we will make the following assumption. 

Hypothesis 2.1. {Vt,F ,v,T) is an invertible probability dynamical sys- 
tem, a : O — > M"' is a v- centered square integrable function and a is a fixed 
real number satisfying < a < 1. 

The invertibility hypothesis is not restrictive since any dynamical system 
admits an invertible extension (its natural extension). Moreover, let us recall 
that any stationary sequence of centered and square integrable random vari- 
ables admits a representation of the form (Y^ = a o T^)k with ($7, F, z/, T) and 
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Fig. 1. 



a as in Hypothesis 2.1. We denote by U'{Vl,W^) the set of measurable func- 
tions f-.n^W^ such that dv{uj) < +oo. For any / E LP(0,M'^), 
we define := {J^\f{u;)\P^di^{uj)y/P. We denote by L°°{n,R'^) the set 
of measurable functions / : ^ R'^ which are i/-almost surely bounded by 
some constant and, for such a function, we denote by ||/||oo the following 
real number: 

ll/IU := inf{Af > : i^{{uj G Q : l/HU > M}) = 0}. 

We consider a system of particles moving independently in W^^^ between 
two d-dimensional horizontal plates M'^ x {0} and M'^ x {e} separated by 
some small distance e > 0. We suppose that these particles move with speed 
i(a(u;),±l) parametrized by ui. In our model, the speed and the parameter 
to only change when the particle hits one of the plates; a particle incoming 
to the upper plate with the speed ^{a{u>), 1) will outgo: 

(a) either (with probability 1 — a) with the speed i(a(T(a;)), — 1); 

(b) or (with probability a) with the speed ^{a{u>'), —1), where u>' is given 
by a random variable (independent of iv) with distribution z/. 

We make analogous assumptions for reflections off the lower plate (replacing 
i by — i and — ^ by ^). We are interested in the behavior of such a model 
when e goes to zero. See Figure 1. 

Let us study the evolution of a single particle moving in this system. Let 
us write 6 = 1 if, at time 0, the particle is pointing upward and 5 = —1 
if, at time 0, the particle is pointing downward. Then, the speed of the 
particle between the nth and the (n + l)st collision off one of the plates is 
i(a(X„),(5(-l)-), where (X„) ^ is a Markov chain such that the conditional 
law of Xn+i with respect to {Xq, . . . , X„) is (1 — a)57^(x„) + ck'^- Let us notice 
that the measure i' is an invariant probability measure for this Markov chain. 
More precisely, we define {Xn)nez as follows. 

Notation 2.1. We consider the probability space (j7,JF, i>) with := 
0^, ^ the product cr-algebra and 9 the unique probability measure defined 
on Q such that we have 

E^[f{Xn)]= I fdv 

Jn 
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and 

Ei>[/(X„+i)|X„,X„_i,...] = (l-a)/(T(X„)) + a / fdu, 

Jn 

with Xn : — > given by Xn{{u>m)m&) '■= ^n- We define the transformation 

r on by f{{u;n)n&) = {iOn+l)n&1- 

Since v is T-invariant, the existence of is a consequence of a result of 
lonescu Tulcea (cf. [21], [32], page 154). With this notation, if the particle is 
at time at the position {x,ez) (with x G M*^ and z E [0; 1]) with the speed 
^(a(Xo), —1) parametrized by Xq, then its horizontal position at time s > 
will be given by 

i~{s,x,z,-):=x + e-z-a{Xo) + £ ^ a{Xk) + e- z\a{X\^^^i^2^^_^^j^i), 

k=i I e J 

where {u} is the fractional part of u. For symmetry reasons, if the particle 
was at time at the position {x,ez) (with x S M'^ and z S [0; 1]) with the 
speed i(a(Xo),l) parametrized by Xq, then its horizontal position at time 
s>0 is given by ^+(s,x,z,-) := ^7(s,x, 1 - z, •)• 

2.2. Results. In [6] Boatto and Golse have studied the following quanti- 
ties: 

F^^{s,x,z,uj) =Ec,[(l){Cf{s,x,z,-))\Xo=uj] 

(their f^^^ corresponds to our F^^ with —a instead of a). More precisely, they 
establish the following result in the situation when the dynamical system 
is given by the algebraic automorphism Tq of the two-dimensional torus 
= ^ given by the matrix (^ ^. We recall that Tq preserves the Haar 
measure vq on T^. 

Theorem 2.2.1 ([6]). Let us suppose that {9.,T,u,T) = {T"^ ,B{T'^),vq, 
Tq). Let (j) be a smooth bounded function with bounded derivatives up to 
order 4. Let a be a u-centered function belonging to H^{T^,M.'^) with (3 > 1 
such that the matrix D (a) :=^f^^^{l — a)^^^'E^[a(S)aoT^] is invertible. Then, 
for any real number tQ> 0, we have 

sup sup sup sup \F^,{s,x, z,uj) — E[(j){x + Bs)]\ = 0{e), 

se[0;to]a;eR'*t<Jen2G[0;l] 

where {Bs)si=u is a d-dimensional Brownian motion with zero mean and with 
covariance matrix D{a). 



RATE OF CONVERGENCE FOR STATIONARY PROCESSES 11 



We will show how this result is related to the central limit theorem for 
(a(Xfc))fc and give some extensions of it. Indeed, for any real number f3 > 1, 
the Sobolev space iJ^(T2,M"') is contained in L°°(T2,M"') and we have 



Remark 2.2.2. Under Hypothesis 2.1, if a is in L°°{Q,W) andcj):] 
is a Lipschitz continuous function, then we have 



sup 

s>0,xeM'',2e[0;i; 



fe=0 



< 



eL0||a||oo(4 + 2^/(l 



a) 



Proof. We have 
For any A; > 1, we have 



+ 1 e 



< 4eL^||a||io 



F;;4(fe+i)e2,x,i,^) = E«'(i-«)'"' E «/o,...^M, 

i=0 «o>0;/i>l,...,/j>l: «o+---+«j=fc 



with 
and 

a«o,...,ii('^) := 



k \ 

m=0 / 



n Jn 



X + e 



V m=0 



+E E 

i=l \ mi=l / 



dl'{uJi) • • • dv(ijjj), 



from which we deduce that we have 

\F"J{k + \)e\x, l,uj) - F-{{k + l)e\x, I, J) 



<2eL<^||a||oo E(^o + l)(l 
/o>0 



a 



Moreover, we have 



f: 



+ 1 e^x,l. 



:E,: 



(t){x + e ^ a{Xk) 



□ 
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Proposition 2.2.3. Let us suppose that Hypothesis 2.1 is satisfied. For 
any integer k>0, we have Ei^[a{Xo)^a{Xk)] = {l — a)^E^[a(g>aoT^]. More- 
over, the following limit exists: 



D(a):= lim Ec 



1 



n-1 



k=0 



and satisfies 



D{a) = - (4^^^y[a ®aoT% 



Proof. This is a consequence of the fact that we have E£;[a(X;i.)|-'^o] 
{l-a)''a{T^{Xo)). □ 

Here, we prove the two following results: 



Theorem 2.2.4 (Rate of convergence in the CLT in the sense of Prokhorov). 



Under Hypothesis 2.1, ifa:i}^'K'^ belongs to L^(0, 



am 



to lL<^/2J+i, 



then {-^ X]fc=o ^i-^k))n converges in distribution to a centered Gaussian ran- 
dom variable with covariance matrix D{a) = J2k^zi^ ~ a)l'^lE^[a (S> ao T'^]. 
Moreover there exists a real number A> such that, for any integer n>l, 
we have 



n-l 



^[i^*[^J2''(^k)],MiO,D{a))\< 



A 



k=0 



'n 

1 v^n-l 



where i^*{-^J2^=o'^i-^k)) denotes the distribution of ^-j^Y^^=QO.{^k) with 
respect to v . 



Theorem 2.2.5 (Rate of convergence in the CLT in the sense of the 
Kantorovich metric). Under Hypothesis 2.1, if a-.Q^W^ is a v-centered 
function belonging to L°°{Q,W^) , then there exists a constant B > such 
that, for any Lipschitz continuous function : R'^ — > M, we have 



(2) 



n-l 



fc=0 



E[0(iV)] 



B , 
<—L4 
n 



where N is a d-dimensional centered Gaussian random variable, centered 
with covariance matrix D{a) = X]fcez(l ~ a)''^'Ej^[a a o T^]. 

Corollary 2.2.6. Under Hypothesis 2.1, if a : il, ^ M.'^ is a v-centered 
function belonging to L°^(r2,]R'^), then, for any Lipschitz continuous function 
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sup 

s>0,x 



<pix + e ^ a{Xk)\ -E[^{x + V^Bi) 

. \ k=0 /. 



where Bi is a d- dimensional Gaussian random variable, centered with co- 
variance matrix D(a). 

Proof. Let a real number s > be given. If s < e"^, then we have 
(j)ix + e a{Xk) 



fc=0 



and 



\&[(I)[x + ^sBi)]-(I){x)\<L^^s\\Bi\\li. 
On the other hand, if s > e^, according to Theorem 2.2.5 [apphed to n = 
[p-J + 1 and to the Lipschitz continuous function z ^ (f){x + ze^^ [p-J + 1)], 
we have 



Ej> (t)\x + e a{Xk) 

. \ k=0 

Moreover, we have 



E 



x + e\ 



+ lBi 



< BUe. 



E 



x + e\ 



+ 1^1 



E[(^(x + ^fsBi 



2^s 



\B 



lIlLi 



^ I^<P2\\Bi\\l^- 



□ 



2.3. Martingale method. We recall that the Markov operator Qa,o asso- 
ciated to is defined by Qafl{f){uj) := E[/(X„+i)|X„ = uj], for any 
/ G L^(0,]R). It is given by the following formula: 

(3) Qa,o(/)(^) = (1 - a)/ o T{uj) + aE,[/]. 

Using a method introduced by Gordin [17], we get 

Proposition 2.3.1. Let us suppose that Hypothesis 2.1 is satisfied. Let 
pG [l,+oo] and a function f G L/{Q,C) u-centered. Then, there exists a 
decomposition of f of the following form: 



f{Xo)=gf{Xo,X_i) + hfiXo) - hf{X_i) 



v-a.s.. 
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with E£>[5/(Xo,X„i)|X_.i] ={) andhe LP{Q,C). 

Moreover, ifp>2, then we have Ej>[(g/(Xo, =D{f). 

Proof. Let p and / be as in the hypothesis of the proposition. Let us no- 
tice that Qafi acts continuously on L^(r2,C) and that we have {Qa,o)^if) = 
(1 — a)""/ o T", for any integer n > 0. Using the fact that o T = X^+i 
(for any k G Z), we can check directly that the functions gf{XQ,X^i) = 

En>0(Qa,0)"(/)(^0) - E,n>l(Qa,0)"(/)(^-l) ^nd hf = - Em>l (Qa,o)" (/) 

are suitable. Let us prove the second point. By definition, we have 
D{f) := lim Ec' ' ' 

n^+oo 

Moreover, we have 

n-l /n-l \ 

Y.f{X,)=( 9fiXk,X,_^) \ + hfiXn-i) - hfiX.i). 

k=0 \k=0 ) 

We conclude by noticing that we have E;>[(/j(Xfc, = if 

The sequence of random variables (X]fc=o5'/("'^fc'^fc-i))" ^ martingale 
and /ij(Xo) — = hfoXo — hfo XqoT~^ is a coboundary in {Cl,D,T). 

This result (Proposition 2.3.1) ensures that, if a is in L^(r2,M'^), then 
the sequence of random variables (a(Xfc))fc satisfies a central limit theorem. 
Here, we are interested in a more quantitative question: the rate of con- 
vergence in the CLT. To this end, we will use more sophisticated methods 
(perturbation of operators. Theorem 1.1). 

Corollary 2.3.2. Let f be a i^- centered function belonging to L^(r2,IR) 
such that D{f) = 0. Then we have / = 0, v-almost surely. 

Proof. Let such a function / be given. According to the two previous 
results, there exists a function h G L^(r2,M) such that 

f{Xi) = h{Xi)-h{Xo), 

Let us show that h is almost surely constant. Let us suppose that there 
exists two disjoint measurable subsets A and B oi^ such that u[h^ A) > 
and v{h £ B) > 0. Then, there exist 0)1,0)2 G ^ such that we have 

Xi{ui)=Xi{u2)=u;, h{uj)eA, h{Xo{ui))eA, /i(Xo(a)2)) £ S 

and 

/(Xi(cD,)) = h{Xi{Cj,)) - hiX^iCji)), 
and therefore /i(Xo(o)i)) = h{uj) — f{ijj) = /i(Xo(o)2)). □ 
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2.4. Proof of Theorem 2.2.4. The idea of the method we present here is 
due to Nagaev [29, 30]. It has been used by many authors (cf., e.g., [19] and 
[20]). From formula (3), we get, for any integer n > 0, 

Qloif){^) = (1 - ")"/ ° + (1 - (1 - ")")E.[/]. 

We recall that we have the following relation: 

Qlo{f){Lo)=E,[f{Xn)\Xo=u;]. 

We will see that the good properties of the Markov operator Qa,o enable us 
to use the method used in particular in [20]. 

Notation 2.2. We denote by (•,•) the usual scalar product in M"^. If 
{B, II • lie) is a complex Banach space, we will use the following notation: 

1. We denote by B' its topological dual (i.e., the set of continuous linear 
maps from B in C). We endow this set of the norm || • ||g/ given by 
\\A\\b> :=sup||j||g=i |^(/)|. 

2. For any AG B' and any / in B, we will use the notation 

{A,f%:=Aif). 

3. For any A G B', any g G B, we denote by g i^^: A the continuous linear 
endomorphism of B defined by 

(5®, A){f) ■.= {A,f),g. 

4. We denote by Cb the set of continuous linear endomorphisms of B. We en- 
dow this set with the norm || • ||£g given by ||-P||£b := sup||j||g=| ||P(/)||e. 

Let us consider the Banach space B :=L°°(r2,C) endowed with the norm 
II ■ IIb = II ■ ||l°°- For any t G W^, we denote by Qa,t the linear operator on 
L^fo, C) (for any q G [1; +oo]) defined by 

QaAf) :=Qa,o(e*<*''^">/(-)). 
With this definition, we have 

(Q„,j(/))(XO=E,[e^<*'"(^"+i))/(X„+i)|X„]. 
The introduction of these operators is motivated by: 

Remark 2.4.1. Under Hypothesis 2.1, for any t G R"' and any integer 
re > 1, we have 

E,[e*<*'i:::o-(^^))|x_i] = (Q,,i)"(i)(^-i)- 
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Proposition 2.4.2. Let us suppose that Hypothesis 2.1 is satisfied. 
Let m>l be an integer. 7/ 0:0 — > M'^ is in L™'(0,M'^), then the function 
Qq,^. : M"^ — > £g is on and, for any integer k = 1, . . . ,m and any 
(ji, ■ • ■ Jk) e {1, . . . we have 

-^i^-TT^Q-^tif) := i^Qa,tM) ■ ■ ■ a,, (•)/(•)) 

(where Oj is the jth coordinate of a). 

We recall that we have Qa,o{f) = + (1 — a)f o T. This can be 

rewritten as follows: 

QaAf) = (1 ®* ^)/ + (1 - o T - E,[/]). 



Theorem 2.4.3 (Perturbation theorem, see [20]). Let us suppose that 
Hypothesis 2.1 is satisfied. Let m > 1 be an integer. We suppose that a : O — > M*^ 
is a V- centered function belonging to L™'(r2,M'^). Then, there exist a neigh- 
borhood Uq of in M.'^ and three nonnegative numbers ci,r]i, r/2 and four 
functions Xa,. G C™(C/o,C), Va,. G C"^{Uo,B), G C"'{Uo,B') and iV„,. G 
C"^{Uo,Cb) such that: 

1. (Initial values) Xao = I, Va o = 1, fao = ctnd Na o(/) = (1 — «)(/ o T — 
E.[/]). 

ox 

2. (Initial derivatives) For any i = 1, . . . ,d, g^;' |t=o = 0; if m>2, then we 

have Hesst A^^f |t=o = -D{a), with D{a) := J2k&z{^ - a)''''lEy [a a o T''] . 

3. For any t in Uq, we have: 

(a) (Decomposition of the operator) For any integer n>l, {Qa,t)"' = 

(b) (Dominating eigenvalue) Qa,tVa,t = Xa,tVa,t, {Qa,t)*^a,t = Xa,tV>a,t 
and {(pa,t,Va,t)* = 1- 

(c) \Xa,t\ > 1 - ?7l. 

(d) For any integer n>l, we have 



max max 

fc=0,...,mii,...,ij.G{l,...,(i} 



d'' 



dti-^ ■ ■ - dt. 



<ci(l-r/i-r/2)". 



Proof. This result is a d-dimensional version of Theorem III-8 of [20], 
page 18. Its proof leads to the implicit function theorem and is exactly the 
same as the proof of Theorem III-8 of [20] . □ 
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2.4.1. Regular case. In our proof, the following theorem plays the same 
role as the Berry-Esseen lemma in the proof of the rate of convergence in 
the one-dimensional central limit theorem (see Theorem B of [20], page 12). 



Proposition 2.4.4 ([44]). Let Q be some nondegenerate d- dimensional 
normal distribution. There exist two real numbers cq > and T > such 
that, for any real number T > and for any Borel probability measure P 
admitting moments of order [|J +1, we have 

U{P,Q) 



< Co 



1 + r 



+ 



\t\oc<T 



Ld/2J+1 

E 

fc=0 {ii 



E 

,ife}e{i,...,d}'= 



Qk 



dti^ ■ ■ ■ dti^ 



{(pp-(PQ){t) 



1/2. 



dt 



We will prove the following and conclude according to Proposition 2.4.4. 

Proposition 2.4.5. Under the hypotheses of Theorem 2.2.4, if D{a) is 
invertible, then there exists a real number (3 > such that, for any integer 
k = 0, . . . , + 1 and any zi, . . . , G {1, . . . , d}, we have 



|iU</3v^ 



gk 



dti^ ■ ■ - dt 



-(E^[e^<*'(^/^)Sr=o''^(^'))] - e-(^/2)<*'^(")*>) 



dt 



1/2 



01^ 

/n 



Proof. The following formula will be useful in the following. 

Lemma 2.4.6. Let k be a positive integer. Let b be a complex-valued 
function C'' -continuous defined on some open subset U ofW^. Letn>l be 
an integer. Let us consider the function u : [/ — > C given by u{t) := [b{-^))"^. 

Then, for any ii, . . . , G {1, . . . , d}, we have 

Qk 

-uit) 



dti^ ■ ■ - dt 



{^1 



n(n ■ 



I) ■ ■ ■ in - m + l)[b[ —= 
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where Qk is the set of partitions A = {Ai, . . . , Am} o/{l, . . . , k} in nonempty 
subsets Ai = \ . ■ ■ il'^^A^}- 

Let C2 > and /? > be two real numbers such that the closed ball 
S|.|^(0,/3) is contained in Uq and such that for any t G 5|.|^(0,/3), we have 
\KA < e-^2(t,t> and e-{i/2)(t,£'(a)t> < ^-C2{t,t) _ ^^r^YAs is possible because D{a) 
is invertible and because we have Hesst Aa,t|t=o = —D{a)?[ In the following, 
n will be any integer and t S M*^ any vector satisfying n > 2 and |t|oo < f3y/n. 
For such a couple we have G Uq. Therefore, we have 

E-[e*(*.(Vv^)Ero'"(^')>] 

= (^,(Qa,Vv^)"l)* 

1. We start by giving an estimation when /c = 0. We have 

IE-[e»(*.(i/v^)Erro''^(^'))] - e-(i/2)(t,D(a)t> 

+ (-^a,t/v^) 



n \Aoo 



n 



+ ci(l -r/i -r/2) 

Therefore, we have 
( I |E,[e*<*'(Vv^) ErJo' «(^'))] - e-(^/2)(t,D{a)t> ,2 ^ A ^ o 1 Y 

2. Let A; be an integer satisfying 1 < < J + 1 and (ii, . . . , i^) G {1, . . . , d}* 
According to Theorem 2.4.3, we have 



c?t,', • • • 9t 



• • • dti^ 



1 
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/ , p-{<^2n){t,t)^ I I Qk 
+ of(i + |i|^;„)-''°'r"U"'^"""'^'" 



((^a,t/v^) 



• • • dti^^ 

+of(i+ie')^''°''""V'''""'""^''' 



since {y,{v^^t/V^ ^* '^a,t/V^W* " 1 = O(^) and gf^^ .g^^^ ((A^.t/VTl)") 
+ |i|^)e~^'^^/^^^*'*^). We will estimate the following quantity: 

^'^ .((X , , )n) -(l/2)(t,D(a)t) 



In the following 6: i?|.|^(0,/?) ^ C will be a function cL'^/^J+i S|.|^(0,/3) 
such that ^(0) = and Hess6(0) = -D{a) and \b{t)\ < e-^2(t,t> [^g take 
b{t) := Aa^t and b{t) := e~(-^/^)^*'^(")*^]. According to Lemma 2.4.6, we have 
9^ / / / t \ \ "\ 

XI 9n,m(^,fe)(i)> 

^={Ai,...,A„}gQfe 



• • • dti^ V V V Vn 
with 



9n,m{A,b){t) := n{n - 1) ■ ■ ■ {n - m + l)[ b[ 



/n 



. ^ \dti , , ■ ■ ■ dti , J \ Jn 

For any ^ = {Ai, . . . , A^} £ Qfc, we denote by mo(^) the number of G ^ 
such that i^Ai = 1. Let us notice that we always have 2m < m^i^A) + h. 
Indeed, we have 

m 

= X > mo(^) + 2(m - mo(^)) = 2m - mo(^). 

i=l 

Therefore, for any A = {Ai, . . . , A^} £ Qk, we have 

/ / Ul \ mo{A)\ 
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= O(n(^/2)(2m-(mo(^)+fe))|^|^'»o(^)g-(c2/2)(i,t>^ 

= 0(|t|S°(-^)e-(^2/2)(M>). 

(a) If ^ = {^1, . . . , Am} G Qk is such that 2m < mo{A) + k, then, for any 
t e i?|.|^(0,/3-y/n), we have 

(b) Now, let us consider a partition ^ = {Ai, . . . , Am} G Qfc such that 
2m = mQ{A) + k. Then, A is made of subsets of {!,..., d} containing at 
most two elements. For such a partition A, for any t G (0, /S-^/n ), we 
have 

\gn,M K-m - 9n,M e-(^/2)<-'^('^)->)(t)| 

= 0(^_L(l + |t|™o(A)+3)g-(c2/2)(M>^_ 



Indeed, we have 



|-(A.,-e-(V2)(-,^W->)f4)=0^" 



|2 

loo 



n / \ n 
3 



and 



.(A„,-e-(i/2)<-.^('^)-))fAUo^"* 



2.4.2. Degenerate case. In this section we suppose that the matrix D{a) 
is degenerate (i.e., noninvertible) . There exists a matrix A G GL{W^) such 
that we have 



A-D{a)-^A = Ji : 



where I is the rank of D{a), Ii is the /-dimensional identity matrix and 
Om,n is the (m, n) -dimensional null matrix. By replacing function a(-) by 
A ■ a(-), we can (and we will) assume that we have D{a) = Ji. According to 
the previous subsection, we have 



(4) n z> J ^ai(Xfc),..., — J^aK^fc) ,AA(0,/0 =0(e), 



where Oj is the ith. coordinate of a. For any i = / + 1, . . . , d, we have D{ai) = 
and therefore, according to Corollary 2.3.2, we have aj = almost surely. 
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2.5. Proof of Theorem 2.2.5. The sequence of random variables {a{Xk))k 
satisfies the hypotheses of Theorem 1.1. Indeed, for any integers a,/3,7 > 
satisfying 1 < a + /3 + 7 < 3, any integers 1< k < k+p < k+p + q < k+p + l, 
and any integers ii,i2,i3, we have 

Ei^[{ai-^{Xk+p))'^ {ai,-,{Xk+p+q))^ {ais{Xk+p+i)y \Xo, Xi, . . . ,Xk] 

=Q%mxk) 

and 

Ep[(aii(Xfc+p))°(aj2(Xfc+p+,))^(ai3(Xfc+p+i))^] =IEi.[V'], 
with ^/' = V'g,/,a,/3,-y,n,i2,i3 :=af^ x (5^^o(«f2 ^ (Q^o^li))- We get 

|Cov(G(Xo, . . .,Xk), {a,,{Xk+p)na,,{Xk+p+,)f{a,,{Xk+p+i)y)\ 
<||Q%(V'-E,[^])|U^||G(Xo,...,Xfc)||ii 
< (1 - - ^miL^ \\G{Xo, . . . 

<2(l-a)P||^|Uoo||G(Xo,...,Xfc)||ii 
<2(l-an|a||S,+'3+^||G(Xo,...,Xfc)||^i 
<2(l-af(l + ||a||^)||G(Xo,...,Xfc)||Li. 

3. Application to the Sinai billiard. The bihiard system considered here 
has been studied in many articles since the fundamental article of Sinai [40] . 
Let us mention [8, 9, 10, 11, 24]. The question of the CLT in this context 
has been studied in many articles [9, 10, 12, 43]. 

Here, we are interested in the question of speed of convergence in the 
CLT. 

A first result has been established in [33] for one-dimensional observables 
(the speed is estimated in the sense of the uniform norm of the difference 
between repartition functions). This result has been extended in [34] (for 
d-dimensional observables, the speed being estimated in the sense of the 
Prokhorov metric). 

The speed obtained in these two papers is in n^^^/^^"*"" for all a > 0. Here, 
we establish a rate of convergence in the CLT in n~^/'^ in the sense of the 
Kantorovich metric. This result is an application of Theorem 1.1. 

3.1. The model. We are interested in the behavior of a point particle 
moving with unit speed in some domain Q of the torus T^, the complement 
of which is a finite union of open sets Oi, . . . ,Oi called obstacles. Each 
obstacle Oi is a strictly convex open set, the boundary of which is and 
the curvature of the boundary is never null. See Figure 2. 

We suppose that the closures of the obstacles are pairwise disjoint. We 
suppose that the point particle moves in Q with unit speed and elastic 
reflection off the obstacles. See Figure 3. 
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3.2. The billiard flow. Let us notice that, when the particle hits an ob- 
stacle, the couple position-speed is ambiguously defined: incoming and out- 
going vectors coexist. To avoid this problem, we decide to take the following 
convention: when a particle hits an obstacle, its position-speed couple cor- 
responds to the outgoing vector. Let us be more precise. For all q in dQ, we 
denote by fi{q) the unit vector normal to dQ in q directed to the interior 
of Q. 

The set of configurations is the set Qi given by 

Qi :={{q,v):q€Q,v€TgQ, \\v\\ = 1 and {q e dQ ^ {n{q) , if) > 0)}. 

We call billiard flow the flow (Yt)t defined on Qi such that, for all t > 0, all 
{q,v) S Qi and all {q' ,v') G Qi, the fact that Yt{q,v) = {q' ,v') means that "if 
a particle is at q with the speed v at time 0, then it will be at q' with speed 
v' at time t." 

This flow preserves the normalized Lebesgue on Qi. 

According to the description of our model, it is natural to study the 
model corresponding to the times when the particle hits an obstacle [cf. the 
system (M, i^, T) below]. 




Fig. 3. 
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Fig. 4. 



3.3. The billiard transformation. Let us consider the billiard system 
(M, I/, T) defined as follows: 

(a) M is the set of configurations of Qi corresponding to the times when 
the particle meets an obstacle, that is, 

M:={{q,v):qedQ,veTgQ,\\v\\=l,{n{q),if)>0}. 

(b) For any i = 1, ...,/, we write /j the length of the boundary dOi of the 
obstacle Oi. We parametrize M by G:M^Ui=i({«} x x [-7r/2; 7r/2]) 
defined by G{q,v) = {i,r,ip) q £ Oi, if r is the curvilinear abscissa of q 
on Oi, and ip is the angular measure taken in [— 7r/2;7r/2] of the angle be- 
tween n{q) and v. 

(c) v is the Borel probability measure on M of the following form: 



where C is some constant. 

(d) T is the transformation of M that, at the configuration {q,v) S M 
of a particle at the time just after a reflection, associates the conflguration 
((?', if)€Mat the time just after the following reflection off dQ (cf. Figure 4). 

(e) We also define the function r : M — > [0; +oo[ where r(g, v) is the time 
to wait for a particle at q with speed v until the next refiection off dQ (cf. 
Figure 4): 




r(g, v) := min{t >0:q + tve dQ}. 
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Fig. 5. 



Let us specify the link between billiard transformation and billiard flow. 
The flow {Yt)t can be viewed as the special flow over the dynamical system 
{M,u,T) associated to the roof-function r. This is very natural: we identify 
{{q,v),s) with {q + sv,v) (cf. Figure 5). 

In the following, we suppose that the billiard system has finite horizon, 
that is, that function r is uniformly bounded. 

In Figure 2 only the second domain corresponds to a billiard system with 
finite horizon. 

3.4. About the regularity of T. Let Rq be the set of configurations cor- 
responding to a vector tangent to an obstacle: 



The study of the billiard is complicated by the discontinuity of T at points 



of T-i(i?o) (cf. Figure 6). 

However, we know that, for any integer > 1, the transformation 
defines a C^-diffeomorphism from M \ Uj^=o ^"^ {^o) onto M \ U^^^q (-^o) • 
Moreover, the sets Uj=o^~"'(-^o) and [Jj=QT^ {Rq) are finite union of C^- 



3.5. Hyperbolic properties of the billiard transformation. For any C - 
curve 7 of M, we define 



using the parametrization of M by the function G previously defined. 



R^:={{q,v)eM:{n{q),v)={)]. 



curves. 
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Proposition 3.5.1. There exist two real numbers Co > and ao G ]0; 1] 
such that, for all x, there exist C^-curves 7'*(x) (stable curve) and j^{x) 
(unstable curve) of M containing x (with positive length for u- almost every 
X ) such that, for any integer n>0, all y,z £ 7*(a;) and all y' ,z' £ l^ix), we 
have 

d{T^{y),T^{z))<Coa^^d{y,z) 

and 

d{T-^{y'),T-'\z')) < Co«^V%',z'). 

3.6. The functional sets 'Hr),m- Because of the discontinuities of T, if 
(f>:M^M. is a Holder continuous function, then (p o T^ is generally not a 
Holder continuous function. This observation leads us to the introduction of 
the sets Tirj^m defined below. These spaces will be such that, if / is r/-Holder 
continuous, then / o is in TCr^^m- 

Let a real number r/ S ]0; 1] be given. For any m, we consider the set Ti.-q^m 
of bounded functions : M ^ C such that the following quantity is finite: 

(jiv,"^) \4>{x) - (t>iy)\ 

■ Sel.x,ylc:x^y (max(d(x,y), . . . , d(r-(x), r-(y))))'' ' 

where Cm is the set of the connected components of M \ 1J^qT~''(/?o) and 
d is the metric defined on each connected component of M by d{{q, v), {q' , if)) = 
a/It - r'l + \(p- (p'\'^ if G{q,v) = {i,r,(f) and G{q',if) = {i,r',(f'). 



r 




Fig. 6. 
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The set W„ m can be understood as the set of functions that are Holder 
continuous in the m future configurations. These classes of functions have 
been introduced in [33]. 

Let us notice that the function r is in Tii^i. 

In the following section, we give a decorrelation result for these classes 
of functions. Before recalling this result, let us make some comments about 
the classes of functions Tirj^m- 

Proposition 3.6.1. Let a real number r/ e]0; 1] and an integer mo > 1 
be given. For any functions <j) o,nd ^ belonging to Ti-r/^mo > '"^^ have: 

1. The functions (j) and ip are uniformly bounded [because the set M \ 
Uj^o-^~"'(-^o) has only a finite number of connected components]. 

2. The function (j) + tp is in Ti-r^^mo '^''^^ have 

3. The product (j) ■ ip is in Ti-ri^mo ^'^^ have 

4. For any integer m>0, (j)o T™ is in 'Hri,mo+m <ind we have 

5. For any integer m>0, the function (p is in 'Hri,m.o+m o-nd we have 

Proof of point 4. Let x and y be two points of M belonging to the 
same connected component of M\[J^^° T~^{Rq). Then T'"(a;) and T^^{y) 

belong to the same connected component of M \ U^o (-^o) and we have 

|<A(T™(x))-c/.(r-(x))| 



< 



Ci"'™"^ max(d(T'"(2;), r™(y)), . . . , ci(r™+"*« (x), T'"+™« (y))"). □ 



3.7. A decorrelation property. The following result has been established 
in [33] (cf. Proposition 1.2 and Corollary B.2 of [33]) with the use of the 
method developed by Young in [43]. 

Proposition 3.7.1. Let a real number r] £]0;1] be given. For any real 
number R> I, there exist two real numbers Crj^R > and 5ri^R G]0; 1[ such 
that, for any integers mi > and m2 > 0, for all (j) G Tirj^m^ and ip G 'Hri,m2! 
for any integer n > 0, we have 

|Cov((/.,V'or")| 

(5) 
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We will not use directly this proposition; we will use a slight modification 
of it: reading the proof of Theorem B.l of [33], we can notice that, in for- 
mula (5), coefficient C^'"^^^ can be replaced by the regularity coefficient of 

on the stable curves: 

rMs))._^^^^ my)-i'{z)\ 

U . .— sup sup — r- 

xeMj;,267=(x) d[y,z)'i 

[by replacing, in the proof of Theorem B.l, the definition of ipk{^) by the 
infimum of ^ o on the stable curve containing x] . 

3.8. Theorem. 

Theorem 3.8.1. Let a real number r] G ]0; 1] and an integer rriQ >1 be 
given. Let / : M ^ R'^ be a bounded function, the coordinates of which are 
in Hr],mo ■ For any k, we write Yfc := / o T'^ and :=Yi + ■ ■ ■ + Yfi. Then, 
the following limit exists: 

T?:= hrn -(E[Sf]). 

n— >+oo n 

// T? = 0, then (Sn) n is bounded in L . 

Otherwise, the sequence of random variables (-^)n>i converges in distri- 
bution to a Gaussian random variable N with null expectation and covari- 
ance matrix and there exists a real number B > such that, for all n> 1 
and all Lipschitz continuous function : — > M, we have 



E 



n 



-mm 



^ BL, 



n 



n 



Proof. For all n, Yi H \-Yn has the same distribution as Xi -\ \-X, 

with Xk '.= f o T~^ . Therefore, it suffices to show that the sequence [Xn 
satisfies the hypotheses of Theorem 1.1. Let us take M := max(l, ||/||oo)- Let 
a, 6, c be integers satisfying a,b,c> 1 and l<a + b + c<3. Let i,j,k,p,q,l 
be integers satisfying l<i<j<k<k + p<k + p-\-q<k-\-p + l. Let 
ii,i2,i3 be integers belonging to {l,...,d}. Let F:W^ x {[-M; M]'^f ^ R 
be a bounded, difFerentiable function, with bounded differential. We have 

|Cov(F(5._i,X,,X,,X,0,(4+i)^4S+.)''(4S+^)')l 

= |Cov(F(/ o T-^ + ••• + / o T"(^-i), / o T-\ / o T-^ / o T-'), 

= |Cov((/>o,V'ooTP+')|, 
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where we define 

and 

iPo ■■= F{f o T^-^+i + ••• + / o T^-i, / o T>'-\ f o /). 
Let us use formula (5) modified with C^'''^*^^ instead of clp'"^^\ 



Lemma 3.8.2. The function -00 

is in T~ij^^2,mo+k—i+i o,nd we have 



i — Oq 



Proof. Let three points x,y,z in M be such that y and z are in 7*(x). 
Then we have 

\f{y)-f{z)\<C^;'""'^CMy,zr/\ 
\f{T^-\y))-f{T^~\z))\<&^'"''\c,a^~yd{y,zf\ 
\f{T^-\y)) - f{T''-'{z))\ < cf>^'^°\Co4-'rd{y,zr/^ 

and 

1/ o r^-*+l(y) + ••• + / o r^-l(y) - (/ o T^-^+l(^) + ... + / o T''-\z))\ 



< |/or'=-+™(y)-/or^-+™(z)| 

m=l 

< ^ cf'"'"\Coa'^-'+"')''d{y,z)''/^ 



m=l 

(fe-i+l)»? 

<C}'''"^0)^.^2 -d{y,zr/\ □ 



Since / is in 7irj,moi the function (po is in 'Hr)/2,mo+i a-^d we have 

||</'0||oo<||/||So+'+^<(l + ||/||L) 

and 

Therefore, according to (5) modified with ci'''^'*^^ instead of ci'''™^^ we have 



|Cov(<Ao,V'oorp+')| 

< C'r;/2,_R(ll</'ol|o< 

<ir(||F|U + ||Z)F|UH,/, 



^ S/2,-Rlll<P0l|oo + <^^o jlll^olloo + ''%/2,R 
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with 

K := C,/2,«(l + 3C}''/^'-"))(l + ml) (l + C^r"^Y^)6~/^-^ 



and 



We have 



APPENDIX 



□ 



Proof of Theorem 1.1. Our proof of Theorem 1.1 follows the same 
scheme as [25]. 

First, since {Xk)k>i is stationary, we notice that, for any integer n > 1, 
we have 



n 



(6) E 

k=l 

Therefore, according to (1), exists and we have 

S2 = E[Xf] + Y.{E[Xi Xk+i] + E[Xfc+i 0Xi]). 

k>l 

Moreover, according to (6), we have 

\nsf] - nS^I^ < 2 J2 k\nXi < 4CM ^ kipk,o- 

k>l k>l 

Hence, if = 0, then the sequence of random variables (5'„)„>i is bounded 
in L2(o,]R'^). 

Let us now suppose that is nonnull. Then, there exists A: S {1, . . . , d} 
and a d-dimensional orthogonal matrix O such that O • • is diagonal. 
Therefore, there exists an invertible matrix A such that 

where denotes the matrix transposed to A and where Jk is the d-dimensional 
diagonal matrix such that the first k diagonal elements are equal to 1 and 
the others to 0. 
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In the following, we will suppose that T? is the d-dimensional identity 
matrix I^i. This is not a restrictive hypothesis: it suffices to replace dhy k 
and {Xn)n by {Xn := {Xn \ ■ ■ ■ ,Xn'^))n where Xn^ is the ith coordinate of 

Xn '■= A ■ Xn [since the random variables {x[''^ -\ \-Xn^)n>i are bounded 

in L^(0,R) for all j = k + 1, . . . ,d and the norms on are equivalent]. 

As in [37] , we will use an inductive proof. The idea is to prove the existence 
of a real number A>1 such that the following property (P„(A)) is satisfied 
for any integer n > 2: 

(P„(A)) : Vfc = 1, . . . - 1, V(/) G Lip(]R'^,IR) 

\E[<j){Sk)] -n^VkN)]] < AL^, 

where is a d-dimensional Gaussian random variable with expectation 
and covariance matrix 1^. Let us define Vn ■= E[5®^] and Vn '■= ^[S'f ^] — 
HS^-i]- We have 

Vn = E[Xf]+ J2 (IE[Ai0A,,+i]+E[Xfc+i0Xi]). 

fc=l,...,n-l 

Hence {vn)n>i converges to = I^. There exists ng > 1 such that for any 
integer n > no, the eigenvalues of Vn are between ^ and |. In the following, 
we will suppose the existence of a sequence (Aj)j>o of independent identically 
distributed Gaussian random variables with expectation and covariance 
matrix Id such that (Aj)j>o is independent of {Xk)k>o- The main part of 
the proof is to establish the following result. 

Proposition A.l. Under the hypotheses of Theorem 1.1, there exist 
a real number K > M and a continuous decreasing function ip: [l,+oo] — > 
]0; +oo[ satisfying limg^+oo ^(s) = such that for any integer n > 9nQ and 
any real number A > M , if we have {Vn{A)), then, for any real number e > 1 
and any Lipschitz continuous function (/) : M'^ — > M, we have 

\E[cP{Sn + sY)] - E[<P{Sno-i + Tno-i,n + eY)] I < K{1 + AiP{e))L^, 

where Y and T„g_i^„ are two v-centered Gaussian random variables indepen- 
dent of {Xk)k>o, with covariance matrices and Vn — Vng-i, respectively. 

Let us show how we can conclude once this result is proved. 

1. Let us write Ai := (i\/9no + maXm=o,...,9no ll'S'mllLi- For any A > Ai, prop- 
erty (VgnoiA)) is satisfied. 

2. Let us show that there exists a real number Aq > M such that, for 
any integer n > 9no and any real number A> Aq, we have (VniA)) =^ 

{-Pn+liA)). 
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Let an integer n > 9nQ and a real number A > M be given such that 
property {Vn{A)) is satisfied. Let (^iR'^ — > M be any Lipschitz continuous 
function. Then, according to Proposition A.l, we have 

|E[<^(5„ + eY)] - E[</>(5„o_i + r„,_i,„ + eY)]\ < K{1 + Ai:{e))L^. 

Since we have 

\E[<i>{Sn + eY)]-n<l^{Sn)]\<L^e¥.[\Y\^] 

and 

|E[(/.(5„„_i + r„„_i,„ + eY)] - E[(A(r„,_i,„)]| 
<L^(eE[|yU]+E[|S„„_i|J), 

we get 

|E[0(5„)] - E[(/)(r„„_i,„)]| < Jfo(l + Ail^{e) + 2e)L^, 
withiro:=i^ + E[|y|oo]+E[|5„o 

— i|oo]- Let us now estimate the fohowing 

quantity: 

|E[0(V^y)]-E[</.(r„,_i,„)]|. 

Let O = On. no be an orthogonal matrix such that OiVn — V^o_i)0 is 
diagonal with nonnegative diagonal coefficients. Let us denote by ^n,no 
the diagonal matrix with nonnegative diagonal coefficients such that 

{l:^n,n,f = 0{Vn-Vn,-l)0-\ 

Let us define M„^„q as follows: 

Mn,no :=0-iA„,„„0. 

Then, we have {Mn^noY = — ^o-i- Let us denote by | • I2 the usual 
Euclidean norm on R*^ and |||j4||| := sup|2|2=i |^2;|2 for any (d, (i)-matrix 
A. Let us recall that if A is a nonnegative symmetric matrix, then |||A||| 
is equal to the maximal eigenvalue of A. We have 

iE[(/>(v^y)] - E[</.(r,„_i,„)]i = \n<t^{v^y) - mn,n,Y)\\ 

< L^¥.[\{^h- Mn,n,)Y\J 
<L^¥.[\{^h-Mn,n,)Y\^] 
<L^\\^h-Mn,nM[\Y\2] 

< L^^\\\nI,-iVn-Vno-i)\m\Y\2] 

< L^.J\\\nId-Vn\\\ + \\\Vn,-l\\\E[\Y\2]. 
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Indeed we have 

\\\^/nId- Mn,no\\\= max |^-A|= max \^/n-^\ 



< , max \/\n-n\ = J \\\nld - {Vn - Vno-i)\l 
neSp{v„-v„Q-i) * * 

where we denote by Sp(74) the set of eigenvalues of the square matrix A. 
Therefore, since the sequence of matrices {Vm — m ■ Id)m is bounded, we 
have 

|E[0(S„)] - E[</.(V^y)]| < K'{1 + Aij{e) + 2e)L^, 

with K' := Kq + E[|y|2] sup^>„Q \/\\mId - Vm\\ + |||Ko-i|||- Let us de- 
note by ea the unique real number ea G [1, +oo[ such that 1 + ^^(sa) = 

EA- 

According to the preceding, for any integer n > 9?io and any real num- 
ber A> M, we have 

{Vn{A)) =^ {Vn+l{3K'EA)). 

Let us show that there exists a real number Aq> M such that, for any 
A> AQ,we have SK'ea < A. The function A>-^ ea^s increasing. If we had 
Ml := sup^ Ea < +00, we would have mi := inf^i ^j^ea) > and therefore, 
for any A, Mi >ea^ A'iP{ea) > Ami, which is impossible. Therefore, we 
have limA—,+oo£A = +oo. Hence, there exists Aq > M such that for any 
^ > ^0, we have 3K'{1 + AiIjIea)) < A and therefore 3K'ea < A. 

Hence, for any real number A > max{AQ,Ai), we have (VgnoiA)) and, for 
any integer n > 9nQ, (Vn{A)) =^ {Vn+i{A)) , from which we deduce Theo- 
rem 1.1. 

Now, we have to prove Proposition A.l. 

Let an integer n > 9no be given. Let (yfc)fc>no be a sequence of indepen- 
dent random variables defined on v) independent of (Xfc)fc>o such that 
Yfc is a Gaussian random variable with expectation and covariance matrix 
ffc. Let y be a Gaussian random variable with expectation and covariance 
matrix Id, defined on {Vt,,T,v), independent of {(Xk)k>nQ-,{Xk)k>o) [this is 
always possible in some extension of {^,J-, v)]. 

Notation A. 2. Let k be an integer such that uq < k < n. We define 
Akif) = E[/(5fc„i + Xk)] - E[f{Sk-i + Yk)]. 

For any function (j) G Lip(M'^,M), for any real number e > 1 and any x G 
W^, we define 



U,k,n,ei^) 

with convention X^iLfc+i Yi = ii k = n. 



\ i=k+l / 
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Let us notice that we have 



E[(l){Sn + eY)]-E 



^k{f4>,k,n,e)- 



We will use Taylor expansions for functions /i:R"'^M. We will use the 
following notation. 

Notation A. 3. Let k be an integer such that k>l. 
If h : Mf^ ^ M is A;-times differentiable, for any x G M'^, we denote by D^h{x) 
the point of Mi^'---''^J' given by 

D^h{x) := ( Kx)) 

We denote by | • |oo the supremum norm on ]Rii'---''^J''° . 

For any A^'^\ . . . , A^^^ in W^, we denote by A^^^ (g) • • • ® A^^^ the point 
of Mi^'---''^} given by 



V=l / ii,...,ik=l, 



[ifA(0 = (4^._^(0)]. 

For any integer j satisfying 1 < j < k, for any A G and any 

B G we denote hy A*B the point of M{i.-..,rf}'^-^ given by 



-( T A 

\ ni,...,nj=l 



'J— J- ' n,...,ifc_j=l,...,d 



For any A:Vl^ Mii.---,'^} ^nd any B:Q,^ ^{'^^■■■A ^ -^q define (when it is 
well defined): 

ll^lloo •- 



Cov{A,B):= Cov{Ai,_i^,Bi,_,^). 

il,...,ife=l 

Let us notice that if j = k, * corresponds to the usual scalar product 
on ]R{1'-''^}'. 

On the other hand, let us notice that the A;-linear form on M'^ associated 
to D^h{x) is 

(Ai, . . . , Afe) ^ D^h{x) *{Ai®---®Ak) 
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and that, for any j = 1, . . . ,k, we have 
{D''h{x)) * (Ai • • • Ak) = {{D^h{x)) * {Ai 
We have Afc(/) = Ai,fc(/) - A2,k(/), with 



)Aj))*{Aj+v 



^Ak). 



and 



A2,k{f) ■= nfiSk^i + n)] - nfiSk-i)] - lE[D^f{Sk-i)] * V,. 



A.l. Estimations for small k. 



Lemma A. 1.1 (Adaptation of Lemma 6 of [37]). Let f : 
function in C^. We have 



be a 



k-l 



|A,,(/)| < d\\\D^f\\^ + pVlloo) + IbC'M'ir + 1) ^(1+^)99^,0 

V p=0 

fc-2 \ 

+ 3CM J2 ^pA- 

p=r+l l=l,...,k-l : {r+l)l<p / 

Proof. Since is a Gaussian random variable independent of Sk~i, 
with expectation and covariance matrix Vk, we have 



|A2,fc(/)| 



E 



f{Sk-i + n) - f{Sk-i) - ^D^fiSk^i) * {Yk ® Yk) 

f\l - tfE[D^f{Sk^i + tYk) * 1^3] dt 
Jo 







{l-t)^E[g{tYk)*Y^']dt 



witlig{u)=E[D^f{Sk^i + u) 
< ^ f\l - tf sup \E[D^f{Su-i + a)]|oo]E[|y®3|oo] dt 

^ Jo naMd 



^3 r\ 

< ^ sup |E[I)3/(Sfc_i + a)]UE[|yf3|^]. 



We have E[|Y'^®^|oo] < Moreover, according to hypothesis (1), we 

have \vk\oo < 2C(M + l)J2pZo'-Pp,o- According to Holder inequality and to 
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the fact that ipp^ < 1, we can show (cf. [37], page 264) that we have < 
(4C7(M + l))3/2'^^J-i(l +p)ippfl. Hence, we have 



|A2,fc(/)|<d^^sup|E[D3/(5fc-i + a)]|oo 
(7) 

k-l 

x(C(M + l)p5](l+p)(^p,o. 

p=0 

Let us now control Ai^fc(/). Since we have Vk = + Ei=/(IE[Xi ® X^] + 

E[Xfc (8)Xj]), we have 

Ai,fc(/) = E[D'f{Sk^,) * Xfc] + i Cov(i?2/(5fc_i), ) 

fe-i 

1=1 

+ E[lD'fiSk-i + ekXk)*Xf% 
where 9k is a random variable with values in [0; 1] . We have 
(8) \\lD'f{Sk-i + ekXk)*Xf'\\^<ld'\\D'f\\^M\ 

On the other hand, according to (1), we have 

|Cov(OV(S»-i),X?")l < X:iCov(OV(S.) - D'f{Si.i).Xf-)\ 

k-l 

<3d^CM\\D^f\\^Y.^pfi- 
p=i 

We have 

k-l 

D^f{Sk-i) = D'fiO) + J2{D^f{S,) - D^fiSi^,)) 

i=l 
k-l 



D'f{0) + J2(D^f{Si^i)*X, 

+ J\^- t)D^'fiS^-l + tXi) * Xf dt 



and so 

k-l 

E[DV(5fe_i) * Xk] - nD^f{Sk-i)] * J2 ^i^^ ® ^k] 



i=l 
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k-l 

= J2Cov{D^f{Si^i),X,0Xk) 

(10) 

fc-i 

+ ^ E[Z)2/(5i-i) - D^f{Sk-i)] * Xk] 
1=1 

k— 1 

+ E Gov ^ (1 - t)D^f{Si^i + tXi) * (Xf 2) dt, X,}j , 

since E[D^f{0) * Xk] = D^f{0) * E[Xk] = 0. According to (1), we have 
k-l 

\nD^f{Si-i) - D^fiSk^i)] * E[Xi Xk]\ 

1=1 

k-l 

(11) <d^D^f\\^M2CMY,{k-i)^k--^,o 

1=1 

k-l 
p=i 



and 



fc-i 



^ Gov f / (1 - t)D^f{S,.i + tX,) * {Xf) dt, X, 
i=i ^-^0 

k-l 

(12) < d'C{2\\D''f\\^M'' + ||Z)Vl|ooM2)(^fc_,,o 

1=1 

k-l 

< 2CdW\D^f\U + \\D^f\\oo)M^ Y Vpfi- 

p=i 

For any integer i = 1, . . . , — 1, we write j = ji := max(0, (r + 2)i — (r + 
According to (1), we have 

\CoviiD''fiSi-i)-D^f{Sj))*Xi,Xk)\ 

i-l 

< Y \Cov{iD^f{S^)-D''fiSm-i))*Xi,Xk)\ 

m=j+l 

<3d^C\\D^f\\ooM\i - j - l)^k-i,o 

and 

mD^fiSi^i) - D^f{Sj)]*E[Xi<S)Xk]\ 
< d^D''f\\^{i - j - l)M2CMipk-i,o. 
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(14) 



Hence, we have 

fc-i 

(13) 

k-l 

< 5d^CM^D^f\\^ir + l)Y.p^p,o. 

p=i 

If (r + 2)i - (r + l)k < 0, then j = and so Cov(i:)2/(5'j), X» ® X^) = 0. 
Hence we have 

k-l 

^\CoyiD^fiSj),Xi^Xk)\ 
1=1 

< E ^|Cov(Z?V(SO-^'/(Si_i),X,®Xfc)| 

i=l,...,k-l: {r+l)k<{r+2)i 1=1 

(r+2)i-{r+l)k 

< d^3CM\\D^f\\^ 

i=l,...,k-l:{r+l)k<{r+2)i 1=1 
k-2 

<3d'CM\\D^f\\^ Y E ^P^r 

p=r+l j=l,...,fe-l: (r+l)i<p ^ 

Lemma A. 1.2 (Adaptation of Lemma 5 of [37]). For any 4> e Lip{M.'^ ,R) , 
for any integer k = uq, . . . ,n and any real number e > 1, the function f(f,^k,n,e 
is C°° and, for any integer i>l, we have 

Ci 

with d := (i*+^3 • 2'''~^ Jj^d |-z|oo|^*^(-z)|oo dz, where h is the density function 
of the Gaussian law with expectation and covariance matrix Id- 

Proof. Let us denote by r„^£2 the positive symmetric matrix such 
that ^ ^2 = K - Vfc + e'^Id- For any x G M"^, we have 



f4>,k,n,e{x) = ^[(I'i^ + ^k+l + --- + Yn + eY)] 

1 



^det{Vn-Vk + e'Id) 



Let an integer i > 1 be given. Let ki, . . . ,ki in {1, . . . ,d} be given. Let us 
denote, for any matrix A, the jth column vector of j4 by [A]j. For any x G M.'^, 
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we have 

(-l)^(Z?7M,n,.(x)),,...,,, 
1 

~^t{Vn-Vk + eH^ 

X <A(n)D^Mr;,t.(n - x)) * (KlJk, • • • KlMdn 

since J^d D^h{z) dz = 0. Let us denote by A_ and A+, respectively, the smah- 
est and biggest eigenvalues of Vn — Vk + sP'Id- Since T^ j^ f,2 is diagonalizable 
in an orthonormal basis, we have 

-J 1/9 

Moreover, for any j = 1, . . . ,d, we have |[r~^^2]j|oo < A_ . Therefore we 
have 

\D' f^^k,nA^)\°o < L<pd^-j^ / \z\oo\D'h{z)\dz. 

Since k>no, we have 

i(n - fc) + < A_ < A+ < f (n - /c) + 

(according to the fact that two invertible symmetric matrices are diagonal 
in a same basis). □ 

Proposition A. 1.3. For any (p G Lip(M'^,M), for any real number e > 
1, we have 

n-[n/3j-l 

E \^k{U,k,n,e)\ < d^Cs + C4)ln(3) 
k=no 

X (m^ + IbC^M^r + 1) ^(1 +p)(/.p,o 

V p>0 

Lp/(r+l)J X 

+ 3CM E E Vp,i]L^. 

p>r+l 1=1 ) 
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Proof. According to Lemmas A. 1.1 and A. 1.2, we have 

n-[n/3j-l 

l^fc(/(/),fc,n,e)| 

fc=no 

<d\C^ + Ci) 

X M3 + 15C2M2(r + 1)^(1 +|5)(^p,o 

m=[n/3j+l + \ p>0 

+ 3CM ^ ^ ^^Al^ 

p>r+l 1=1 J 

A.2. Estimations for big k. 



□ 



Lemma A. 2.1 (Analogous to Lemma 7 of [37]). For any real number 
^> /o?^ o.ny integer n > 9no, if property {Vn{A)) is satisfied, then, for 
any (j) G Lip(R'^,R), for any real number e > 1, for any integer k £ [n — 
LfJ ; ^l; /o'" '^'^2/ integer / G ; A:], /or any integer i>l, we have 

(15) sup |E[«'/,,.„,(5, + < /f. + 

OTi/i i^i := 2^/2(i^ /igd 11)^/1(^)100 + CiVf~^ . 

Proof. Let Ni be a Gaussian random variable with null expectation 
and covariance matrix I ■ 1^. Let ^2 be as in the proof of Lemma A. 1.2. 
First, let us notice that, if {Vn{A)) is satisfied, we have 

|IE[I)7<^,fc,n,e(5, + O)] - E[D'f^^k,nANl + «)]|oo 

- E[<j){Ni + a + T^^k,e^ ■z)]\- \D'h{z)\oo dz 



< 



-j^ AL,\D^h{z)\^dz, 



[n — k + 

since y ^ (j){y + a + F„ ^2z) is Lipschitz continuous with Lipschitz constant 
bounded by L^. On the other hand, we have 

nD'f^^k,n,e{Ni + a)] = D\a^{nf4>,k,n,e{Ni + a)])) 
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where is a Gaussian random variable with null expectation and covariance 
matrix Id- As in the proof of Lemma A. 1.2, we get 



Ci 



yjn-k + l + e^ 
According to estimation (7) and since /c — 1 > we have: 



n 



-La,. 



□ 



Lemma A. 2. 2. Let an integer n > 9no and a real number A > M be 
given. If property (Vn{A)) is satisfied, then, for any (j) G Lip(IR'^,]R), for any 
integer k = n— [^J , . . . ,n, for any real number e > 1, we have 



I A. 



2,k\Jih,k,n,e) 



(16) 



4l(V2^ 

k-l 
p=0 



A 



+ 



(n — A; + e^)'^/^ nj 



iVcM)3/2 



Lemma A. 2. 3. There exists a real number K (only depending on d, C , 
M and r ) such that, for any (j) G Lip(M'^, R), for any integer n > 9no and any 
real number A > M, if property ('P„(yl)) is satisfied, then, for any integer 
k = n— [^J ,n, we have 



\^l,k{f<j),k,n,e)\ <L^K 



(17) 



Aa 



n—k 



+ 



n—k 



(re — fc + e^)'^/^ n 



n~k 



+ 



with am-=l + T}p=iPCp, Pm-=l + Ep=i Cp, Im ■= Ep=|„^/(r+2)2] Cp and 



Proof. To simplify notation, we will write fk instead of f(f,^k,n,e- We 
have 

Ai,fc(/fc) =EpVfc(5fc-i) * Xk] + iE[Z)Vfc(Sfe_i) * (Xf - Vk)] 
+ \nD^fk{Sk-i)*Xf] 



/V - tfD^fj,{Su^^ + tX^) * (Xf 4) dt 
Jo 
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Since we have Vk = ^Xf^]+J2^=l,...,k-lO^[Xk'S)Xi] +E[X,(g)Xk]) and since 
the matrix E[D^/fc(S'fc_i)] is symmetric, we have 

Ai,k{fk)=nD^fk{Sk-i)*Xk]- J2 nD^fk{Sk-i)]*nxk0x^] 

fi=i,...,fe-i 

+ iCov(I)Vfc(Sfc_i),Xf ) 



+ iE[D3/fc(5fc_i)*Xf] 



/ {l-tfD'fk{Sk^i + tXk)*XfUt 
Jo 

For any integer i £ {0, ... ,k — 1}, we have 



: Xi {D^fk{S,)-D^h{S,^^)) 
j=k-i 

k-l 

: D^f,iS,.l)*X, 

j=k-e 

+ i J2 D'fkiSj-i)*Xf 
j=k-e 

k-l 



+ 5 E / {l-tfD''fk{S,., + tX,)*Xfdt. 



Therefore we have 



Ai,fe(/fc)=E 



i f\l-tfD^fk{Sk-i+tXk)*Xfdt 
Jo 



+ lE[D^fk{Sk-i)*Xf^] 
+ E[Z?Vfc(5fc_,_i)*Xfc] 

- nD''fk{Sk-i)]*nxk®Xi] 

i=l 

+ ^Cov{D^fk{Sk-i),Xf^) 

k—l 

+ i Y nD^fk{S,^i)*{Xf^Xk)] 
j=k~e 



42 



F. PENE 



+ 1 E IE [ /\l - tfD^fkiSj^i + tXj) * (Xk 05 Xf) dt 
k-l 

- J2 nD^fkiSk^i)-D^fk{Sj.i)]*E[Xk(^Xj] 
j=k-e 

k-i 

+ J2 Cov{D^fkiSj.i),X,0Xk). 
j=k-e 

In the following, we take £ = [y/n — k\, the integer part of y/n — k. Since 
n > 9no > 9 and ^ < /c < n, we have Q<l<k-l. 

A.2.1. Control o/E[i /q (1 - tfD^f\{Sk-i + tXk)Xf^ dt]. According to 
Lemma A. 1.2, we have 







E 


\-f 




.6 Jo 



-f {l-tfD'MSk-i + tXk)*XfUt 

D JO 



(18) 



<d^-\\D^fk\\ooM' 



<d' 



a 



'n—k 



24 (n- A; + £2)3/2 



A.2.2. Control of IE[D^ fk{Sk-i) * Xf]. We have D^fkiSk-i) 
D^fk{Sk-i-i) + EUi(D''fkiSk-j) - D^fk{Sk-j-i)). According to (1) and 



to Lemma A. 1.2, we have 

■i3 



(19) |Cov(I)Vfc(Sfc-^-i),^r)l <^%^^(C3 + C4)9'Lv^j+i,o 
and 



J2\CoviiD^MSk-,) - D^h{Sk^.j^i)),Xf 



(20) 



3C4M 



^ (n - A; + £2)3/2 
Aa 



< d-'Ls^CiM- 



n—k 



(n-fc + £2)3/2' 

Moreover, according to Lemma A.2.1, we have 



nD^fk{Sk-i)]*nxr]\ 



(21) 



{n — k + £2)3/2 ~^ n} 
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A. 2.3. Control of W.[D^ fk{Sk-i-i) * X^]. According to (1) and to Lem- 
ma A. 1.2, we have 

\E[D^fk{Sk-i-i)*Xk]\ = \CoY{D^fk{Sk-i-i),Xk)\ 

(22) 



A.2.4. Control of E^+i |E[Z)2/fc(Sfc-i)] * nXk-j ®Xk\\. We have 



fc-1 

\nD''fk{Sk^i)]*nxk~j®xu]\ 

k-l 

<(f \nD''fk{Sk^i)]\oo\nxk^j(^xi,]\, 

(23) 



k-l 



6^ 



A.2.5. Control of \CoY{D^fk{Sk-i),Xf^)\. We have 

e 

D^fkiSk-i) = D^fkiSk.e.i) + Y{D^ fk{Sk-j) - ^'/fc(Sfc_,-i)) 

i=i 

= D'^fkiSk^e^i) + Y ^'/fc(5fc~i-i) * Xk^j 
i=i 



e 1 

+ {l-t)D^fk{Sk^j^i + tXk-j)*X^^.dt. 



i=i 

Hence, we have 

Cov(Z)2/fe(Sfc-l),Xf2) 

i 

= CoviD^fkiSk-e-i),Xf) + Y CoviD'fkiSk-i-i) * Xk-j,Xf) 

i=i 

e e 

+ E E CoviiD^fk{Sk-m)-D^fkiSk-m-l))*Xk-j,Xf^) 
j=l m=j+l 
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i ,1 



+ ^ / (1 - t) Cov(Z?Vfc(Sfc-,-i + t^fc-,) * ) di- 

1. First, according to (1) and to Lemma A. 1.2, we have 

|Cov(Z)Vfc(5fe-£-i),^f )l < d'C^£iS=L<^^,+i,o 
(24) ^ + ^ 

2. Control of ^^=1 Qo^ip^ h{Sk-i-x) * Xk-j,Xf). 

(a) For any integer j = satisfying ^"^"^ < i < ^, we have 

\/n — /c < (r + 2)j and, according to (1) and to Lemma A. 1.2, we have 



Hence 



e 

^ \CoY{D^fk{Sk-e-i) * Xk-j,X^ 

j=[v/^^/(r+2)J+l 

(25) 



32C(C3 + Q)M 

(b) For any integer j = 1, . . . ,£ satisfying j < ^'^"^ , according to (1) 
and to Lemma A. 1.2, we have 

|Cov(Z)3/fc(5fc-£-i),^fc-,®^f )l<d3^%^^0mi-,,j 

and 

|Cov(2?3^fc(5fc_,_i),Xfc„,) *E[Xf ]| < ^3 g^'(g3 +C'4) 

Moreover, according to (1) and to Lemma A. 2.1, since k — i — 1 > ^, we 
have 

|E[D'A(St-»-i)l • E[^i-, » ^ril 



from which we get 

L\/^/(r+2)J 



5] |C0Vp3^fc(5fc_,_i)*Xfc_,-,Xf ) 
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(26) < (f2CM^{C^ + C4) 2 ^'f> 

3. Control otE5.iEi.^.+iCov((D3A(St-m)-0='A(Sl-».-l))»^t-,j,^f: 

(a) For any integers j and m satisfying 1 < j < m < i and m < (r + 2) j , 
according to (1) and to Lemma A. 1.2, we have 

\CoviiD^fkiSk-m)-D^fkiSk-m-l))*Xk-j,X^'')\ 

Hence, we have 

£ min(^,(r+2)j) 

^ J2 \CoviiD^fk{Sk-m)-D^MSk-m-l))*Xk-j,Xf)\ 
j=l m=j+l 

(27) 

(n- A; + £2)3/2 ^<^- 

(b) We will use the following formula: Cov(yl * B,C) = Cov{A,B 
C) - Cov{A,B) * E[C] + E[A] * E[B (C - E[C])]. For any integers j 
and m satisfying 1 < j < (r + 2)j + 1 < m < ^, according to (1) and to 
Lemma A. 1.2, we have 



^ 3C4M ^ 

and 

|Cov(Z)Vfc(5fc-m) - DVfe(Sfc-™-l),Xfc_,-) * E[Xf ] 



3C4M3 



'(n + £2)3/2 ^'^'^—i.O' 
from which we get 



£ 

E \CoviD^fkiSk.,n) - D^fk{Sk-m-l),Xk-j )| 

i=l (r+2)i+l<m<£ 

+ E E |cov(d3/,(5,_^) 

i=l (r+2)j+l<m<£ 
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(28) 



6CCaM^ 



(r + l)(n-A; + e2)3/2 

Indeed, we have 

e 

E E 



D''fk{Sk^m^i),Xk^j)*E[Xl 



a 



j=l {r+2)j+l<m<i 



e lp/(r+l)i 

< E E '^p.i^vrf 

p=r+2 j=l 



and 



E E 'Pm-jfi < 

i=l (r+2)j+l<m<^ 



a 



'n—k 



r + 1 



Finally, according to (1) and to Lemma A. 2.1 (since k — i — 1> |^), we 
have 



E 



5] E[Z;Vfc(^fe-m,)-^'/fe(5fc_ 
m=min(£,(r+2)j)+l 



E[Xfc_,55(Xf _E[Xf])] 



(29) 



4. We have 



)-I)Vfe(5fc-£-l)] 

*E[Xfc_,®(Xf -E[Xf])]| 

L£/(r+2)J 



L^/(r+2)J 



< dMCK-,M 



A 



(n — A; + e^)-^/^ n/ 

H ]-Lrl 



+ -]L^2CMipj^o 



(n - A; + £2)3/2 



n 



^ /\l - t)|Cov(Z)Vfc(Sfc-,--i + tX,„,) * Xf ^,Xf )| 
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(30) <rf^E [\l-tf-^^^^±^^L,^,odt 

.7 = 1 ■ ■ 



(n-fe + £2)3/2 
A; + £2)3/2 «V^^ 



^ ,4 3C(C4 + C5)M2 



A.2.6. Contra/ o/ Ej=i E[L'3/fe(5fe_j_i) * (X®^ ® X^)]. Let us notice 
that we have 

i e 
J2 nD^fkiSk-j-i) * . (^Xk]=Y^ Cov{D^fk{Sk-j-i) * Xf^.,Xk). 
i=i i=i 

We control this quantity as we did for 

e 

i=i 

in the previous section (we get analogous estimations). 

A.2.7. Control of Ej=i/o(l-i)'lE[Z)Vfc(5fc-,-i + tXfc_,)*(Xf ^.^X^)] dt. 
We control quantity Y!j=i /o (1 Cov(Z)V-fc(5fc_,_i + tXfe_,) ., X^) dt 
as we did for E5=i/o (1 " t) Coy {D^ h{Sk-j-i + tX^-,) * Xf^.,Xf)dt. 
We obtain estimations analogous to (30). 

A.2.8. Control of E,'=iIEp'/fc(Sfc-i) - I?Vfc(5fe-,-i)] * E[Xfe_, ® Xk]. 
For any integer j = 1, . . . , we have 

fk{Sk-i) - fk{Sk-j-i) 

3 

= E (-D fk{Sk-m) - D fk{Sk-m-l)) 

= X] fk{Sk~m-l) * Xk-m 

■m=l ^ 

+ - t)I)V(5fe~m-l + tXk-n.) * Xfl^dt^ . 

Moreover we have 

D^fk{Sk-m-l) = fk{Sk-l-l) + E {D^ fk{Sk-p) - fk{Sk-p-l)). 

p=m+l 
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Y,nD''fk{Sk-l)-D^h{Sk-j-i)]*nXk-j®Xk] 



i j 



J2T.^ f {'i--t)D^fk{Sk^m-l+tXk-m)*Xf^,^dt 
■ 1 1 



j=l m=l 



: E[Xk_j ® Xk] 



I j 

nD^fk{Sk-e-i) * Xk^rn] * nxk-j ® Xk] 

j=l m=l 

+EE E nD\MSk-p)-D^fkiSk-p-i))*Xk. 

j=l m=l p=m+l 

*E[Xk^j®Xk]. 



1. According to (1) and to Lemma A. 1.2, we have 



j=l m=l 



E E ^ ^ (1 - t)D^fk{Sk-m-l + tXk-m) * Xf^^dt 

*E[Xk-j(^Xk 



(31) 



<^^EE 



(n-A; + £2)3/2 



<d 



4 2CM^Ci 

(n - fc + £2)3/2 



2. We have Ej=i Ei=i E[DVfc(5fc-^-i) * ^fc-m] * ^Xk-j ® Xk] = + As 
with 

{£>(r+2)m}) 

j=l m=l 

{£<(r+2)m}- 

jr = l m=l 
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According to (1) and to Lemma A. 1.2, we have 

j=l m=l 

< d'CM^ uJ2 ^^ E + 1 - "^)mi-m,o 

rn=l 

(32) <d'CM^^^±^L, ± p^,,o 

p=[(r+l)VS^/(r+2)] 

(33) <d^CM^^^^L,,^. 

[We use the fact that if i> {r + 2)m, then we have m < ^'^"^^ and £ + 

1 _ m > > + l)m and m < 

On the other hand, we have 

I . I / ^3V- V- 2C7M2C3 

l^2| <a 2^ 2^ I £2'^i.ol{^<('-+2)m}^0 

j=l m=l 

,s 2CM2C3 



o 2CM2C3 

3. We have 



(35) 



and 



e j {r+2)j 

\n{D'fk{Sk-p)-D'fkiSk-p^^))*Xk^rn] 

j=l m=l p=m+l 

<EE E (,_,^,.)3/. w^^.o 

j=lm=lp=m+l ^ ' 

2(i4CC4M3 
- 7 — 9T^('^ + 2)a,A73fc-t^(/. 



E |E[(i?-Vfe(5fc_p) - i)Vfc(5fc_p_i)) * 

j=lm=l p={r+2)j+l 
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*E[Xk-j'^Xk]\ 

(36) <'^'EE E 7rrTT7w^'^^"° 

i=lm=lp=(r+2)i+l ft + ej/ 

(n- A: + £2)3/2 (r + ipE^^ '^^'"^^ 

^ ,4 3C4M3 1 

- (n - A; + £2)3/2 + 1)2 "v/^^^^</" 

A. 2.9. Control of Y!j=iGoY{D'^fk{Sk-j-i),Xk-j®Xk). For any integer 
J = 1, we have 

D'^ fk{Sk-j-l) = D'^ fk{Sk-e-l) + E {^"^ fk{Sk-m) — D"^ fk{Sk-m-l)) 

m=j+l 

and, for any integer m = j + 1, . . . ,i: 

D'fk{Sk-rn)-D'fk{Sk —rn—l ) 

= D^fk{Sk-m-l) * Xk-m + /V - t)D^f{Sk-m-l + tXk-m) * Xf^^dt 

Jo 

e 

= D^fk{Sk-£-l) * Xk-m+ E i^^fk{Sk-p)-D^fk{Sk-p-l))*Xk-m 

p=m+l 

+ f\l- t)D^f{Sk-m-l + tXk-m) * X^^^dt. 

Jo 

Therefore, we have 

e 

J2 Cow{D^fk{Sk-j-i),Xk-, ® Xk) 

£ 

= E Cow{D^fk{Sk-£-i),Xk-j ® Xk) 

+ E E Gov / {l-t)D^fk{Sk-m-l+tXk-m)Xf^^dt, 
j=l m=j+l ^-^0 

Xk-j eg) Xfc 

+ E E Cov(DVfc(5fc-£-l)*^fc-m,^fc-i«'^fc)l 
j=l m=j+l 
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{e>(r+2)m} 

j=l m=j+l 
I I min{€,{r+2)m) 

j=l m=j+l p=m+l 

+ E E Cov((I)Vfc(5fc-(r.+2)m-l)-^Vfe(5fc„£-l)) 
j=l m=j+l 

* Xk-m,Xk-j Xk)l{£>(^r+2)m+l}- 

Let us control each term of the right-hand member of this identity. 

1. [Control of ^$=1 Cov(L>2/fc(5-t-^-i), ^fc-i ® ^fc)-] 
According to (1) and to Lemma A. 1.2, we have 

e 

^ \CoviD^fkiSk^i^i),Xk^,0Xk)\ 

i=L€/(r+2)J+l 

< Yl i\Cov{D''fk{Sk-e-i)*Xk_,,Xk)\ 

j=[e/{r+2)\+l 

and 

L£/(r-+2)J 

^ |CovpVfc(5fc-^-i),^fc-j®^fc)l 

W{r+2)} ^^^^ ^ 

Let us notice that if j < Lt^^J j then we have (r + 2)j < i and so j = 

^ < ^ and > ^- L42J +1 > fi^V^ > There- 

fore, we have 

e 

J2\Cov{D^fk{Sk-e-i),Xk-j(^Xk)\ 
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^ ,2 ^C{C2 + C^)M ^ 



\/n — k + e"^ ^ 
2. [Control of E5=i E^=i+i Cov(/o^(l -t)I)Vfe(>5fc-,„-i + *Xf ^ dt, 

Let j and m be two integers satisfying 1 < j < j + 1 < m < £. 
(a) If ?n < (r + 2)j, then we have 

Gov - t)D^h{Sk-m-l + tXk-m) * {Xflm ® ^k-j) dt, Xfc) 

t)D^fk{Sk-m-l + tXk-,n) * X®!^ dt 

Jo 



E 



■■E[X, 



k-j 



I Xu 



^ , 3{C, + C,)M^ 
(b) If ?n > (r + 2)j, then we have 

C0V^\l - 

Therefore we have 

7 = lm=7 + l ^ ° ^ 



(38) < 



{n- A: + £2)3/2 



E f'v'p,o+ E E v^pj)^</ 

\ p=l p=l j=l / 



- („_A: + e2)3/2^</" 

for some Kq only depending on d, C, C4, C5, M and r. 

3. We have 

e I 

J2 E |Cov(Z?3/fc(5fe„,„i)*Xfe_^,X,._,®Xfe)l 

j=l m=j+l 



<d 



,3 



n — k + 



/ LVn-fcJ 

E ^"^P.o 

V=rVn^/(»-+2)21 
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(39) ^ 

p=|-(r+l)V7i^/(r+2)2] i=l / 

n — k + e"^ 

for some only depending on (i, C, M and r. Indeed, let j and m be 
two integers satisfying 1 < j < j + 1 < < ^. 

(a) If £ < (r + 2)m and m < (r + 2)j, then we have 

72 rC ~r £ 



and 

(b) If < (r + 2)m and m > (r + 2)j, then we have 



, 2(C3 + C4) M 
n — + e 
4. We have 

{e>{r+2)m} 

j=lm=j+l 

r+2)m} 

j=lm=j+l 

+ 11 51 Cov(D3/fc(5fc_^_i),Xfc_„®Xfc_j 

{£>(r+2)m} 

i=l ■m=j+l 

e e 

-Y H nD^fk{Sk-i^l)*Xk-rn]*nXk-j^Xk]l{e>^ 
i=l m=j+l 

(a) We have 

j=l m=j+l 

+ . 2)3/2+^)^0 
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Since we have 

I (r+2)j £ (r+2)j 

\n^k-m®Xk-j®Xk]\oo<Y. E 2cmVj,o 

j=lm=j+l j=lm=j+l 

I 

<2CAf2^(r + l)jVj-o 
i=i 

and 

i i it 

j=l m=(r+2)j+l i=l m=(r+2)j+l 

£ b/(^+i)J 
p=i j=i 



we get 

j=l m=j+l 

(40) 

73/1/ I ^\m\Jt'i■J^ ( ^^\fn-k , PVn-k 



\{n — k + e^'y/'^ n J 
(b) For the two other terms, we write 

C3 + C4 

|Cov(i:>^/fc(5'fe_^_l),Xfc_.„ < (fC —-^L^ipi+i_rn,m 

and 

Ti Ki ~\~ S 

according to (1) and to Lemma A. 1.2. Let us notice that, if £ > (r + 2)m, 
^ and so£ + l-m>^^:±^ 



then we have m < ^" „^ and so £ + 1 — m > _ > (r + l)m and 



'TT- < ■ Hence, we have 



'/2'£+l-m,ml{£>(r+2)m} 

j=lm=j+l 

l-m,ml{£>(r+2)m} 

m=2 
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< > max 7 

Therefore, we have 

X! XI |Cov(L»^/fc(5fe_£_i),Xfc_^0Xfc_j ®Xfc)|l|^>(^+2)m} 
j=lm=j+l 

(41) 

In the same way, we get 

i £ 

{£>(r+2)m} 

(42) ^'^'"^^^'^^ 

n — K + 

5. We have 

I e min{e,{r+2)m) 

EE E \CoviiD^fk{Sk-p)-D^MSk-p-i)) 

j=l m=j+l p=m+l 

(43) 

-(n-fe+V^4 E E P E 

- (7i-fc + £2)3/2 

for some -fCa and A'3 only depending on d, C, C4, M and r. Indeed, let j, 
m and p be three integers satisfying l<j<j + l<m<m + l<p<i. 

(a) If we have p < {r + 2)m and m < (r + 2)j, then we have 



(44) 
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and 

\E[{D^fk{Sk-p) - D^h{Sk-p-i)) * Xk-m] * nXk-j Xk]\ 

(b) On the other hand, ii p <{r + 2)m and m > (r + 2)j + 1, then we 
have 



\Cov{{D-'fkiSk^p) - D-'fk{Sk^p^i)) * Xk^m,Xk-j Xk)\ 

(n - A; + £2)3/2' 
We conclude with the use of the fohowing formulas: 



e (r+2)i (r+2)m £ 
j=lm=j+lp=m+l j=l 

and 

e e (r+2)m I I 

E E E V'm-ij=E E (r- + iWm-j- J- 

j=l m=(r+2)j+l p=m+l j=l m=(r+2)j+l 

£ Lp/(r+i)J 

<(^ + 2)E E V^v^r 
p=i j=i 

6. Now we control 

^ CoviiD^MSk.ir+2)m^l)-D^fk{Sk^i^l)) 
j=lm=j+l 

* Xk-m,Xk-j (8) ^/c)l{£>(r+2)m+l}- 

If ^ > (r + 2)m + 1, then we have 

r+2)m-l) ~ D fk{Sk-e-l) 

e 

= iD^fkiSk-p)-D^fkiSk^p^i)). 

p=(r+2)m.+l 

We will use the following formula: 

Cov(yl *B,C(g)D) = Cov{A, B®C®D) 

- E[A * * E[C (g) D] + E[^] * E[5 (g) C L>] . 
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(a) If £ > p > (r + 2)m + 1, then we have 

^ 3C4M ^ 

and 

\E[{D^fk{Sk-p) - D'USk-p-i)) * Xk-m] * nXk-j ® Xk]\ 

The sum of these quantities over {j^m,p) satisfying l<j<j + l<m 
and (r + 2)m + 1 < p < £ is less than 



(n-fc + e2)3/2- 



4\AI^J [VH^j Lp/{''-+i)J 

\ p=l p=l j=l 



Thus we have 

I I 

Y Y \Cov{D^fkiSk-(r+2)m-l) 
j=lm=j+l 

(46) - -D^/fc(S'fc_^_i),Xfc„.m (g) Xfc-j 18) -'^fc)|l{£>(r+2)m+l} 

-\n-k + e^)y^ 



and 



^ ^ |E[(Z?Vfe(Sfc_(.+2)™-i) 

j=lm=j+l 



(47) 



* E[Xk~j Xk]\l{e>(^r+2)m.+l} 



- (n-fc + e2)3/2 ? 



(b) Let us now control the following quantity: 

e e 

^ E[D''h{S,_(^r+2)rn-l)-D^fkiSk-i-l)] 
j=l m=j+l 

*E[Xfc„m (g) Xk^j Xk]l{i>{r+2)rn+l}- 
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If £ > (r + 2)m + 1 and m > (r + 2)j + 1, we have 

|E[0'A(S/l-(r+2)m-i)l • E|(Xi_„ - E[Xi_„J) » Xt-j ® Xill 

If £ > (r + 2)m + 1 and m < (r + 2)j, we have 

|E[OVt(Si,-(,+2)m-l)l • nXt-,n ^l-j 8) (Xt - E[Xt])]| 

Therefore, we have 

e e 

{^>(r+2)m+l} 

j=l m=j+l 

^ ^^^6 f 7 , ^ 2^3/2 + - 

V(n — fc + e^)'^/"' n 

\ p=i j=i p=i / 

- ^ ^^(n-fc + . 2)3/2 + ^^j^^- 

In the same way, we get 

e I 

j=l m=j+l 

(49) 



(48) 



, (n — + £2)3/2 n 
This completes the proof of Lemma A. 2. 3. □ 

A. 3. End of the proof of Proposition A.l. Let an integer n > 9no and 
a real number A> M he given. Let us suppose that property {Vn{A)) is 
satisfied. Let a real number e > 1 be given. 

Let us recall that we have 

n 

E[(t>{Sn + SY)] - E[^{Sno~l + Tn,-l,n + sY)] = E ^k{U,k,n,e), 

k=no 
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with TnQ-i^n ■= J27=no Then, according to Proposition A. 1.3, we have 

n-[n/3j-l 

\^kiU,k,n,e)\ < K'-IL^, 

k=no 

where K'^ only depends on {d,C, M,r, {ipp^i)p^i). 

On the other hand, according to Lemma A. 2. 2, we have 



(50) 



/ 1 



k=n— [n/3j 



Moreover, according to Lemma A. 2. 3, for any integer k = n— [^J , . . . , n, we 
have 



\^l,k{f<j,,k,n,e)\ < L^K 



(n- A; + £2)3/2 



n 



n — + 



Ad 



+ 



n~k 



with am = 1 + Ep=i Kp, /3m = 1 + Ep=i Cp, 7m = Ep=|m/(r.+2)21 and 5m = 

T.p=\rn/{r+2)\+i ^ ^1*^ Cp := P maxj=o,..., Lp/(r+i)J V'pj- We control indepen- 
dently each sum of these terms over A; S {n — [|^J , . . . , n}. 

1. Control of the first term: 
" Aa 



E 



n—k 



k=n— [n/3j 



(n- A; + £2)3/2 



(51) 



L«/3J 1 


= ^E- 




1=0 


<A\ 


((^ 


<A\ 


((^ 



(/ + £2)3/2 
(/ + £2)3/2 



p?aS(^^^ 

E . ^ 



Kp 



pCp 



2. Control of the second term: 



(52) 



E 



k=n— [n/3j 



n 



<i + ECp- 

p>i 
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3. Control of the third term: 



fc=n-Ln/3j «=0 v=\^^l|{r+2Y'\ 

Lv^Jp2(r+2)4 

p=0 i=p2 



Therefore, we have 

'''' .Xj^^' "^-'^^^ 

4. Control of the fourth term: 

n as: L"/3J -I +00 

y A^VTT-^ ^ ^ y ^_ y Cp 

fe=n-[n/3j i=0 ^ ' p=|^^/(r+2)J+l 

p2(r+2)2 



<A5]2(v/p2(^ + 2)2 + e2_V?TT 

Hence, we have 



(54) E -7=^^^E 



/ ,V-2(r + 2)2(1 +p2)^^ 



,=„_L„/3j Vn-fc + e2 Vp2 + e2 

5. Control of the fifth term: 

n L"-/3J 

E '^LV^J+i-o= E '''Lv^J+I.O 

k=n-[n/3\ 1=0 

<^#{/:LV/J=p}(/.p+i,o, 

p>0 

from which we get 

n 

(55) E <^Lv^J+i,o^E(2^'+l)^p+i,o- 

fc=n-[n/3j p>0 
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